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Abstract 

Nitrogen vacancy (NV) centers in diamond have distinct promise as solid-state qubits. This is because of their large 
dipole moment, convenient level structure and very long room-temperature coherence times. In general, a combination 
of ion irradiation and subsequent annealing is used to create the centers, however for the rigorous demands of quantum 
computing all processes need to be optimized, and decoherence due to the residual damage caused by the implantation 
process itself must be mitigated. To that end we have studied photoluminescence (PL) from NV~, NV" and GRl 
centers formed by ion implantation of 2MeV He ions over a wide range of fluences. The sample was annealed at 600° C 
to minimize residual vacancy diffusion, allowing for the concurrent analysis of PL from NV centers and irradiation 
induced vacancies (GRl). We find non-monotic PL intensities with increasing ion fluence, monotonic increasing PL in 
NV°/NV- and GR1/(NV° + NV^) ratios, and increasing inhomogeneous broadening of the zero-phonon lines with 
increasing ion fluence. All these results shed important light on the optimal formation conditions for NV qubits. We 
apply our findings to an off-resonant photonic quantum memory scheme using vibronic sidebands. 
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1. Introduction 

The creation of single photoluminescence cen- 
ters in diamond has promising applications in the 
development of hybrid optical solid state quan- 
tum devices. Diamond exhibits a large inventory of 
optically active centers, related to the vibrational 
and electronic states of impurities and defects in 
the crystal lattice [1]. The prime candidate for 
many such applications is the negatively charged 
nitrogen- vacancy (NV^) center. It has proven to 
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be a reliable source for single photons [5] suitable 
for quantum key distribution [5], and when used 
as a spin qubit has demonstrated coherent oscilla- 
tions [3], qubit tomography Stark shift control 
of the optical transition and coherent population 
trapping [6j. Coupling of NV~ centers to nearby 
nuclear [S] and electron [718] spins has also been 
demonstrated. Furthermore, technology to sculpt 
optically important nano and micro-structures in 
diamond is progressing rapidly [5]. Based on such 
promising results, diamond seems a strong candi- 
date platform for quantum information processing, 
and several different schemes have been proposed, 
including resonant dipolc-dipole coupling, repeat- 
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until-success and brokered graph states [TU]- A 
recent review of this field can be found in Ref . [llj . 

The schemes mentioned above were all based on 
the formation and control of isolated qubits, how- 
ever there are many protocols which benefit from 
ensembles of optically active centers. In ensemble 
protocols, the emphasis is usually on using the sys- 
tem as a medium to act on individual photons or 
photonic qubits, rather than focusing on the cen- 
ters themselves as the qubits. Here we are explicitly 
thinking about protocols such as the weak nonlinear 
coupling schemes |12| . or photon storage schemes 
[13ll4j . Again, NV~ appears to be an extremely 
suitable scheme for such applications, and electro- 
magnetically induced transparency has been demon- 
strated at radio-frequency [H] and optical transi- 
tions [16j in ensembles. Our principal concern in the 
present work is the optimization of the formation 
methods of ensembles for these applications, with 
special focus on storage of light outlined in Ref. [14] . 

Beyond the obvious attraction of ensembles in in- 
creasing the density of particles over the single par- 
ticle case, ensembles of emitters have other advan- 
tages. The electromagnetic interaction strength is 
strongly enhanced [17] , which increases the fidelity 
of state transfer and storage of quantum informa- 
tion |18j . Of great pertinence to quantum memory 
schemes is that the distributed storage of quantum 
information provides more temporal stability. More- 
over, ensembles account for the modal structure of 
interacting light and allow for the storage of the tem- 
poral shape of qubit photons [M] . 

The realization of scalable devices for the above- 
mentioned applications requires adequate control of 
the formation processes of active NV centers in di- 
amond crystals. As the final applications of the NV 
centers differ from the previous motivations (which 
were usually into the fundamental properties of the 
centers), so too are the formation requirements. In 
particular, for quantum device applications, we re- 
quire the ability to engineer spatially defined ensem- 
bles, perhaps in some patterned array, with mini- 
mal residual inhomogeneity, and compatibility with 
sculpted optical structures. These are demanding re- 
quirements indeed, and force a critical reappraisal 
of the existing methodologies. 

At least two alternative routes are available for 
the creation of NV centers in single crystal diamond 
(a third, where NV centers are created during dia- 
mond growth is more pertinent to diamond formed 
by chemical vapor deposition, e.g. [H]). Firstly radi- 
ation damage from various sources (electrons, neu- 



trons, ions, etc.) can be used to create vacancies in 
nitrogen-rich (i.e. type lb) diamond, which combine 
to form NV post-anneal. Secondly, direct implanta- 
tion of nitrogen in the purest (i.e. type Ila) crys- 
tals, where the NV center is formed by the vacancies 
due to the N implantation. While the latter strategy 
represents a suitable method to fabricate low den- 
sity NV ensembles or isolated qubits in diamond, 
the first strategy provides a ready method to quickly 
and efficiently produce high density ensembles. Be- 
cause of our focus on ensemble protocols, this is the 
approach that we have chosen in this study. Here, 
we present an extensive study on the effects of MeV 
He+ ion implantation in type lb diamond crystals, 
as a candidate technique to create high-density NV 
ensembles in diamond. Irradiation with high energy 
ions creates vacancies as the ion looses energy to the 
lattice. Such processes can be easily and accurately 
modeled using the standard TRIM package [20] . Be- 
cause the energy loss of the ions is non-uniform with 
distance into the material, with most damage cre- 
ated at the end of range of the ion, this perforce 
leads to a highly non-uniform depth profile of NV 
centers formed post-anneal. 

Although NV ensembles in diamond are promis- 
ing for optical quantum information processing 
(QIP), there are several material-related issues that 
need to be fully understood if ion irradiation is to be 
considered for the creation of NV center ensembles: 

NV center density. A high conversion efficiency 
from nitrogen and vacancy defects to optically ac- 
tive NV centers is crucial to reach high NV densi- 
ties, increasing the probability of an interaction, but 
also minimizing the number of unconverted N which 
will contribute to decoherence without contributing 
to the interaction. An active NV center can be cre- 
ated in a diamond crystal containing nitrogen and 
vacancy defects with thermal annealing at temper- 
ature > 600°C. At such temperatures [H], the va- 
cancies start migrating to the nearest substitutional 
nitrogen atoms, where their aggregation is energeti- 
cally favorable 22] ■ The conversion efficiency is lim- 
ited due to competitive processes such as the forma- 
tion of other defects and vacancy-interstitial recom- 
bination. 

Temporal stability of the NV center. The 

NV defect can exist in two charge-states (NV~ and 
NV°) that have been correlated with A=638 nm 
and A=575 nm luminescence emissions, respectively 
[1]. Charge transfer [23l24l25l26j mechanisms lead 
to the conversion between NV~ and NV° centers 
(photo-ionization) . This process appears experimen- 
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Fig. 1. Typical emission scheme for a color center in diamond 
for short wavelength excitation and at cryogenic tempera- 
tures. After excitation, fast non-radiative decay transfers the 
population to the vibrational ground state of the electronic 
excited state. From here, radiative decay occurs both into the 
zero phonon line (ZPL) and the phonon sidebands (PSBs). 

tally as photo-bleaching. The equilibrium between 
NV~ and NV" concentrations is determined by the 
presence of nitrogen donors in the lattice |27I28| . 
Only the NV~ center has an easily accessible A- 
shaped energy desired for many quantum computa- 
tion applications, and a high NV^ / NV° ratio is 
desirable. 

Inhomogeneous spectral broadening. The 

inhomogeneous broadening of the zero phonon line 
(ZPL) NV~ emission is usually reported as 750 
GHz |T], which significantly exceeds the homoge- 
neous linewidth corresponding to a transition life- 
time ~12 ns. Such large broadening is commonly 
attributed to variations in strain and electric fields 
within the diamond crystal, due both to impurities 
and structural defects [TT] , although recent work in 
low N diamond where NV concentrations less than 
the canonical 750 GHz have been observed seem 
to indicate that this linewidth can be altered by 
suitable preparation. Inhomogeneous broadening 
can lead to dephasing or undesired transitions in 
resonant QIP schemes, while off-resonant [T^ or 
frequency-selective [6] schemes are less vulnerable 
to it. 

The above mentioned issues are related with the 
fabrication of NV centers in diamond. They might 
impose limits on the suitability of NV centers for 
certain QIP schemes. In the present work, we com- 
pare our findings on NV center ensembles with 
the requirements for a photonic quantum memory 




Fig. 2. Off-resonant storage scheme of a broadband pho- 
ton q into a local phonon mode of a defect center, for ex- 
ample, the NV~ center. A denotes the detuning of the 
transition from resonance, C is the classical control pulse. 
The actual storage is achieved in an ensemble of ab- 
sorbers. 

scheme (Qmem) that uses broadband (i.e. sub- 
picosecond) pulses and an off-resonant transition 
|14j . often called Raman quantum memory |29] . 
This scheme utilizes a strong classical control pulse 
to absorb the qubit photon into a storage state in 
the material, which is a collective state of the en- 
semble. Unlike the typical QIP schemes that use the 
spin-state of the NV center, here, the first excited 
phonon state has been chosen as the storage state 
(see figure[2]). Most optically active defect centers in 
diamond show phonon sidebands (PSBs) [SDIj and 
PSB storage as introduced here could in principle be 
realised in many defect centers. The thermal popu- 
lation of this state is 10~^° at T ^ 4 K, making state 
preparation, e.g. by optical pumping, unnecessary. 
The read-out of the stored photon can be achieved 
using another strong classical control pulse. 

Raman-like transitions are typically insensitive to 
inhomogeneous broadening of the upper level, pro- 
vided that the energy separation between ground 
and storage state is well-defined throughout the en- 
semble. Certainly this is the case for the phonon 
sidebands considered here. However, with increas- 
ing inhomogeneity, a larger detuning is necessary 
and a stronger control field is necessary to reach a 
sufficient strength of the two-photon transition. On 
the other hand, a large detuning allows for the off- 
resonant use of truly broadband pulses. For the NV~ 
center, the storage state splitting is = 15.3 THz, 
allowing for a bandwidth oi 5 ^ 3 THz, i.e. sub- 
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picosecond photon storage. Due to the high refrac- 
tive index of diamond, the storage is optimum for ir- 
radiation at the Brewster angle, to avoid reflection- 
loss of the qubit photon. 

Our scheme is not affected by the strong phonon 
sideband emission typical for the NV center, as the 
transition is truly off-resonant and the excited state 
is never populated. This is obvious at storage, as the 
transition from state | 2) to | 3) occurs stimulated 
by the strong control field and thus only into the 
one desired PSB. At readout, emission of a single 
NV center would occur into all phonon sidebands, 
which would lead to a non-deterministic frequency 
of the recreated photon. While this affects Qmcms 
based on single NV centers (and additional mea- 
sures have to be taken such as suppression of PSB 
emission by microcavities) , this problem is automat- 
ically avoided by an ensemble quantum memory. 
Due to collective enhancement of ensemble emission, 
the qubit photon is emitted at the ZPL frequency 
only. For any other transitions (PSB transitions), 
no phase-matching |31| can be achieved, effectively 
suppressing the emission into these modes. 

Taking into account the practical considerations 
for light storage in a realistic sample of NV dia- 
mond, we find that at a flucncc F = 1 x 10^^ ions 
cm~^, a A = 648 nm photon with picosecond band- 
width could be stored using a 16 fiW ultrafast con- 
trol laser at 80 MHz repetition rate (transition oc- 
curs at AA = 10 nm detuning from the ZPL reso- 
nance). With these parameters, the storage and re- 
trieval fidelities would range in the percent region, 
and a proof of principle experiment could be estab- 
lished. Higher Qmem fidelities require higher optical 
depths, which could be achieved by different implan- 
tation techniques, e.g. multiple energy ion implanta- 
tion or high energy electron implantation. It has to 
be stressed that, although we calculate this scheme 
for the NV~ center only, it can in principle be ap- 
plied to any optically active defect in diamond with 
a A-shaped energy level structure, for example the 
GRl defect, which also has an optically accessible 
phonon-sideband. This implies that these phononic 
schemes for light storage arc considerably more flex- 
ible than the more conventional approach using ex- 
isting atomic or molecular levels. The current sam- 
ple therefore has been annealed at a comparatively 
low temperature (see next chapter) , to allow for the 
PL-analysis of both NV and GRl. 



2. Experimental Setup 

2.1. Implantation and Annealing 

The employed sample was an artiflcial diamond 
produced by Sumitomo. The crystal was cut and pol- 
ished from a large single crystal which was synthe- 
sized under high pressure and temperature (HPHT). 
All sides of the sample have [100] face orientation, 
and the crystal is classified as type lb; the nitrogen 
concentration is ~ 100 ppm, as reported by the man- 
ufacturers. All irradiated regions were in the central 
[100] growth sector of the sample, where negligible 
fluctuations in the nitrogen concentration can be 
expected, as verified by cathodoluminescence map- 
ping. 




'■■ " Depth {^im) 

Fig. 3. (a) Transmission optical microscope image of the 
implanted sample; the 80x80 fim^ implanted regions are la- 
beled with the relevant ion fluences (ions cm~^). Regions 
implanted at higher fluences are clearly visible because the 
damaged layer is opaque, while regions implanted at fluences 
below 10^^ ions cm~^ are not visible optically, (b) TRIM 
Monte Carlo simulation of the damage density profile in- 
duced in diamond by 2 MeV He"*" ions, scaled for an implan- 
tation density of 1 X 10^*^ ions cm^^; the damage density 
is concentrated at the end of the range, namely ~ 3.5 fira 
below the sample surface. 

The sample was implanted with 2 MeV He+ ions 
on the MP2 microbeam line of the 5U NEC Pel- 
letron accelerator at The University of Melbourne. 
Regions of 80 x 80 /im^ were implanted in the central 
100 growth sector of the sample at fluences ranging 
from 1 X 10^^ ions cm~^ to 2x10^^ ions cm~^, us- 
ing a raster scanning ion beam which was focused 
to a micrometer-sized spot, providing homogeneous 
ion fluence. The ion beam current ranged from ~5 
pA to ^^0.5 nA for different implantations, with im- 
plantation times < 10 — 30 min. The ion ffuence 
was measured by monitoring the Rutherford back- 
scattered ions, after having coated the sample with 
carbon in order to avoid charging; afterwards, the 
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conductive coating was removed by chemically etch- 
ing. The sample was tilted by an arbitrary angle in 
order to avoid channeling effects. 

FigurelHK shows a transmission microscope image 
of the sample after ion implantation; the regions im- 
planted at higher fluences are clearly visible because 
the high damage density makes the material opaque, 
while at fluences below 1x10^^ ions cm~^, the im- 
planted areas are invisible to optical microscopy. 
Figure[3b shows the vacancy density profile of 2 MeV 
He"*" ions in diamond, calculated with "Transport 
of ions in matter" (TRIM) Monte Carlo simulation 
code [23. The simulation takes into account both 
direct collisions and secondary recoil events, and the 
value of the atom displacement energy was set to 50 
eV. Ion-induced damage causes a progressive amor- 
phization of the crystalline structure with increas- 
ing ion fluences, which is localized mainly at the end 
of range of implanted ions (namely ~ 3.5 fim below 
the sample surface, as shown in Fig.tSb), where most 
of the nuclear collisions occur. 

The evolution of the amorphized structure upon 
thermal annealing depends critically on the dam- 
age density. In regions where the vacancy density is 
above a critical threshold, the diamond crystalline 
structure is permanently converted to a sp^-bonded 
phase. When the damage density is below the critical 
threshold, thermal annealing has the effect of con- 
verting the amorphized structure back to the crys- 
talline diamond phase, although residual point de- 
fects can form in the crystal. During the thermal an- 
nealing, vacancies (interstitials) start becoming mo- 
bile in the crystal lattice at temperatures of ^ 600° C 
[21j : in this process, they can cither rccombinc with 
interstitial atoms, or form a vast range of lumines- 
cent centers, both isolated and paired with native 
impurities (N, B, Si, Ni, etc.) [J. 

The sample was annealed for 1 hour in forming 
gas (4% hydrogen in argon) atmosphere. The sam- 
ple was annealed at the relatively low temperature 
of 600°C, which represents the onset for vacancy mi- 
gration in diamond. This low annealing temperature 
induces only a small mobility of the isolated vacan- 
cies, avoiding a strong reduction in their concentra- 
tion, and only a small fraction of them migrate to 
substitutional nitrogen atoms where they combine 
to form NV centers. This allows for the concurrent 
presence of both NV~ , NV" and GRl centers in the 
material, whose photoluminescence analysis is the 
aim of this work. Using these TRIM calculations, 
the He ion fluence has been converted into the aver- 
age vacancy density of the 3-/.(m-thick cap layer be- 



fore annealing. After the annealing less than 0.25% 
of the vacancies is expected to anneal out [32j , thus 
making our estimation of the damage-induced va- 
cancies reasonably accurate also after the annealing 
processing. The cap layer is defined as the first 3 
/im of the sample, as shown in figure [3b. While this 
calculation is a good estimate at low ion fluences, it 
has to be noted that it is approximative only at high 
fluences due to non-linear material processes scaling 
(e.g. graphitization) . 

2.2. Optical Experiments 

PL experiments were carried out using a Rcnishaw 
Raman apparatus with a Leica DMLM microscope. 
The sample was cooled to T ~ 4 K in an Oxford In- 
struments liquid helium flow cryostat, which could 
be positioned with micrometer precision. A reflect- 
ing microscope objective was employed for both ex- 
citation and signal collection, minimizing the effects 
of chromatic aberration. The lateral resolution of 
the system was ~ 2 /xm and therefore much smaller 
than the size of the implanted regions. The photo- 
luminescence was induced by Aoxc = 532 nm light 
from a Suwtech single longitude mode laser. 




550 700 850 1000 

Wavelength (nm) 

Fig. 4. PL spectra for various irradiation fluences. The ZPLs 
of the color centers have been indicated by arrows, their FSB 
are to the right. The left end of the spectra, containing the 
NV" peak, is magnified by a factor of 5. The vacancy density 
(cm~'^) due to ion implantation in the cap layer (see text) 
is denoted on the right. 

All spectra were acquired with the excitation and 
collection focus on the sample surface. Since the 
depth proflle of the induced damage is non- uniform, 
slightly different focusing conditions might in prin- 
ciple affect the intensity of the PL signal; normal- 
ization to the intensity of the first order diamond 
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Raman line is not possible sinee this is also signif- 
ieantly affeeted by ion-induced damage. Therefore, 
the reproducibility of the focusing conditions were 
tested by repeated measurements at various sample 
positions, from which variations of only 2% were ob- 
served. Our spectra were recorded using the spec- 
trometer with 300 and 1200 grooves/mm diffraction 
gratings. With Aoxc = 532 nm laser excitation, the 
NV° PL emission and first order Raman line (Ra- 
man shift=1332 cm^^, corresponding to ARaman = 
572 nm) partially overlap. The Raman emission has 
been subtracted from the NV° PL signal. While an 
absolute estimate of the color center densities from 
the photoluminescence signal is difficult and error- 
prone, we can derive relative trends of the center 
densities using the relative intensities of PL signals 
acquired at the various implanted regions of the sam- 
ple, which is necessary for optimizing the conditions 
for NV formation. 

3. Results 
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Fig. 5. Intensity of the NV", NV° and GRl zero-phonon— 
line PL as a function of the implanted ion fluence (upper X 
scale) and corresponding vacancy density in the cap layer 
(lower X scale; see text for definition). The NV and GRl 
emission intensities have a non-monotonic behavior, exhibit- 
ing maxima at 3 X 10^^ vacancies cm""' and 6 X 10^" vacan- 
cies cm"'' for NV and GRl, respectively. 



Figure|4]shows PL spectra from regions implanted 
at 14 different ion fluences, ranging from 2 x 10"'^^ 
ions cm~^ to 1 x 10^^ ions cm~^. The NV~ and 
GRl peaks (A=638 nm and A=742 nm, respectively) 
are clearly visible for most regions, although the 
GRl doublet structure is not resolved due to an 
overlap of the two peaks attributed to inhomoge- 
neous broadening. The NV*^ emission is at A=575 
nm and much weaker, but visible for some regions. 
The broad peaks on the red side of each ZPL peak 
arise from vibronic emission (phonon side bands). 

3.1. PL intensities as a funetion of ion fluence 

Fig. [5] shows the intensities of the ZPL emissions 
from NV^ , NV" and GRl centers as a function of the 
implanted ion fluence. In all cases, the PL intensities 
follow a non-monotonic trend. After a steady sub- 
linear increase, the NV^ and NV° emissions reach a 
maximum at fluence F = 1 x 10^^ ions cm"^, then 
steadily decrease until maximum fluence F = 2 x 
10^^ ions cm~^ is reached. In interpreting this trend, 
it is necessary to consider the number of vacancies 
created by implanted He"'' ions at various fluences, 
as compared with the native substitutional nitro- 
gen concentration of ~ 2 x 10^^ cm"'^, correspond- 
ing to a 100 ppm concentration. By integration of 
the vacancy profile reported in Fig. [31d, one can es- 
timate that each 2 MeV IIc+ ion on average creates 



38 vacancies in the diamond crystal. The vacancy 
density profile is strongly non-uniform with depth, 
with ~ 80% of the vacancies being created in the last 
0.5 fim of the total 3.5 fim ion damage range (high- 
lighted in Fig.[3}3). As visible in the transmission mi- 
croscope image in Fig. [3^, regions implanted at flu- 
ences F > 1x10^'' ions cm~^ turn opaque because of 
graphitization in the buried heavily-damaged layer 
3.5 /im below the diamond surface. Therefore, at 
high ion fluences only the ^ 3 /im-thick surface layer 
is responsible for the enhanced PL signal, the end- 
of-range layer being optically opaque. 

At fluence F = 1 x 10^'^ ions cm~^, about 3 x 10^^ 
cm""^ vacancies are created in the surface layer with 
a relatively uniform depth distribution. As this va- 
cancy density exceeds the nitrogen density ('^ 2 x 
10^^ cm"^), the increase of NV PL with increasing 
irradiation fluence stops at F > 1 x 10^^ ions cm~^. 
However, due to the low annealing temperature, not 
all vacancies associate with nitrogen atoms (see be- 
low). 

A conflrmation for this interpretation can be 
found in the increase of the relative intensity of the 
GRl emission with respect to the NV emission, as 
shown in Fig. [6l The GRl emission is unambigu- 
ously attributed to isolated neutral vacancies in the 
diamond structure [T], and such a trend clearly in- 
dicates that isolated vacancies progressively exceed 
the NV centers population. It is worth remarking 
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that GRl emission is measured at low (i.e. < 1 x 10^'^ 
cm~^) ion fluences; this indicates that even where 
the native substitutional nitrogen concentration is 
higher than the vacancy concentration, not all the 
created vacancies get trapped at an N atom during 
thermal annealing. Again, this is due to a limited 
efficiency in the N+ V^NV conversion process [12] , 
as is expectable for the relatively low annealing 
temperature used. 

As shown in Fig.[5l the NV^ and NV'' intensities, 
after having reached a maximum, start decreasing at 
higher ion fluences. While at first a plateau behav- 
ior might be expected after the saturation of native 
nitrogen at F ~ 1 x 10^^ ions cm~^, it must be con- 
sidered that the increasing damage density is detri- 
mental to the NV emission, both in terms of the in- 
creasing optical absorption of the damaged crystal 
and of the competitive formation of other lumines- 
cent centers. Damage induced absorption is qualita- 
tively visible in the optical microscope image (figure 
[3^) from regions implanted at fluences higher than 
F = 1 X 10^^ ions cm~^. 

At high ion fiuences. the strong suppression of 
NV PL detected indicates that graphitization starts 
affecting even the cap layer. Although it is likely that 
the NV density increases in the material, most of the 
NV centers are not optically accessible any more. 

For the annealing process employed here, the op- 
timum damage density for the formation of optically 
accessible NV ensembles in type lb diamond by He 
ion implantation is ~ 3 x 10^^ vacancies cm~'^; such 
a density can be achieved with medium ion implan- 
tation fluences (^^ 1 x 10^'^ cm~^ for 2 McV He+ 
ions). As a higher annealing temperature would in- 
crease the vacancy mobility and therefore revert the 
material damage to some extend, a higher anneal- 
ing temperature is expected to shift the NV peaks 
in figure [5] to the right. Annealing temperatures of 
800° C have been reported to maximize the density 
of optically active NV centers [25] , while the GRl PL 
would be strongly reduced at higher annealing tem- 
peratures [32 • The GRl emission also shows a non- 
monotonic trend, with a maximum at F = 2 x 10^^ 
ions cm~^. Again, the abrupt decrease of the GRl 
intensity at higher fluences can be explained with 
the damage-induced optical absorption. 

Assuming a ^ 0.25% efficiency in the conver- 
sion from nitrogen- vacancy pairs to optically active 
NV centers upon thermal annealing, we can esti- 
mate a density of optically active NV centers of 10^^ 
cm"'^ in the cap layer. This corresponds to an optical 
depth of D = 0.01, while for higher annealing tem- 



peratures D ^ 0.2 could be reached due to a higher 
conversion efficiency, corresponding to a number of 
centers of N-^y- = 5 x 10^ in a 1 fj,m confocal spot. 
Here, the optical depth has been calculated from 

1 dp 1'12U 

D = 

2heQC 7 ' 

where 1^12 denotes the ZPL frequency, 7 the de- 
cay rate of the excited state, a the surface den- 
sity of absorbers, and d the dipole moment (esti- 
mated from the radiated lifetime Ti = 12 ns as c? = 
[Sneohc^ -ff 12 /u^]^^^ ~ 5.5 x lO'^o Cm, where /12 
is the fraction of state | 2) decay that goes into the 
ZPL). 
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Table 1 

Qmem fidelity rj of storage or retrieval for an NV" ensem- 
ble Qmem, with necessary cooperativity parameters C and 
microcavity Q-factors. An implantation density oi F = 10^^ 
ions cm^^ and a N — ► NV conversion efficiency of 5% for 
the cap layer have been assumed. 

In the context of the off-resonant Qmem scheme 
introduced above, a high ensemble density is 
necessary to reach a high fidelity. A generic result 
derived by Gorshkov ^3] covering all three types of 
Qmems (Raman, FIT and CRIB) states that the 
fidelities for storage and retrieval each are given by 

C 



where C ^ D is the cooperativity parameter, and 
D the optical depth. Thus, for NV densities corre- 
sponding to a fiuence of F = 1 x 10^^ ions cm^^, and 
optimised annealing temperature of 800 degrees, a 
Qmem would have an efficiency of 77 17%, which 
is sufficient for a proof of concept. Practically, suffi- 
cient control field energy would be easy to achieve. 
Assuming an 80 MHz repetition rate and a micro- 
PL setup, the transition could be realized for a 
control-field laser average power of 1 mW at a AA = 
10 nm detuning from the ZPL resonance. However, 
for realistic QIP applications, a much higher stor- 
age fidelity is necessary, which for ion implantated 
samples could be reached by a longer propagation 
length, e.g. through a waveguide-structure [3], or 
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Fig. 6. Ratios of photoluminescence intensities. At about 
3 X 10'^'^ vacancies cm~^, the fraction of NV~ decreases due 
to a lack of free negative charge in the material, provided 
by substitutional nitrogen atoms (sec text). 

by an enhanced light-matter interaction due to a 
microcavity. Here, a microcavity of much lower Q- 
factor would suffice, e.g. Q ^ 1000, compared to 
Q-factors several orders of magnitude higher that 
are necessary for high fidelity Qmems based on 
single NV centers. 

3.2. Ratio ofNV- and NV" PL 

Few studies have been conducted on the relative 
NV"/NV° intensity as a function of induced dam- 
age, using radiation sources (neutrons, electrons) 
that do not allow a straightforward estimation of the 
induced damage. In the context of QIP with NV~ 
centers, a very high ratio is desirable, as the occur- 
rence of NV" disturbs optical transitions involving 
NV~ centers due to a spatial variation of the re- 
fractive index and temporal instability. This will be 
explained in more detail in the next section. 

A common characteristic of the spectra reported 
here is the low intensity of the NV" emission with 
respect to the NV~ and GRl emission. This is due 
to the fact that the type lb diamond sample has a 
relatively high concentration of single substitutional 
nitrogen. Nitrogen acts as an electron donor in di- 
amond with an ionization energy of 1.7 eV, while 
the ground state of of an NV^ center lies 2.58 eV 
below the conduction band [34] . The electrons pro- 
vided by nitrogen donors are responsible for the con- 
version of nearby NV" centers to the charged NV~ 
state [2S]; the process can be described by the re- 
action: NV° N ^ NV" + N+. Studies carried 
out on samples with different concentrations of sub- 
stitutional nitrogen proved that at low concentra- 



tion, the NV emission arises mainly from the neutral 
centers, while at high concentration, NV~ lumines- 
cence becomes largely predominant |27I4| . This is 
explained by statistical arguments with a variation 
of the Fermi level as a function of N doping concen- 
tration |35|24|25] , or microscopically with the prox- 
imity of the NV defect to a nitrogen donor |28I26| . 
Our NV~ /NV° ratio of about 10 for low and medium 
implantation densities is consistent with previous re- 
ports [27], measuring an NV~/NV" ratio of 3.3 for 
a CVD grown sample with a medium-high nitrogen 
density of 47 ppm. As the photoluminescence excita- 
tion efficiency most likely varies for NV" and NV~ , 
a constant factor applies to this ratio for different 
excitation powers. 

This, however, docs not affect the relative trend in 
Fig.m which shows the evolution of the NV~/NV° 
ratio with increasing implantation fluence. Although 
the NV^ and NV" emissions follow a similar non- 
monotonic trend with respect to the nitrogen satura- 
tion point at F = 1 X 10^^ ions cm^^^ ^he NV"/NV° 
ratio steadily decreases for increasing fluences, both 
in the low fluence and high fluence regimes. This 
behavior can be explained with the charge transfer 
process that determines the NV° <^ NV^ conver- 
sion. With increasing damage density, more nitro- 
gen atoms pair with induced vacancies to form NV 
centers and fewer nitrogen donors are available to 
turn NV° centers into NV~ centers. This imposes 
a natural limit to the NV~ density achievable with 
high energy He ion implantation for a given sort of 
diamond. As the presence of NV° is unfavorable for 
QIP, no more NV centers must be created than ni- 
trogen atoms are available to provide the negative 
charge for it, i.e. in the case of full conversion effi- 
ciency, the vacancy density created should not ex- 
ceed half the density of nitrogen of the sample, leav- 
ing an electron donating nitrogen atom for each NV 
center formed. However, due to the limited mobility 
of the vacancies at the annealing temperature of our 
work, this saturation effect arises only slowly and 
at a higher vacancy densities of 3 x 10^" vacancies 
cm^'^. 

3.3. Photochromism 

Since the attribution of the 575 nm emission to 
the neutral state of the NV defect by Mita [35] , con- 
clusive work has demonstrated the photochromic 
behavior of the NV center [36124126] . Photoioniza- 
tion has been suggested as an explanation for pho- 
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Fig. 7. Excitation power dependence of the NV PL (normal- 
ized to excitation power) and PL intensity ratio (vacancy 
density 1.5 X 10^* vacancies cm~^). Photoionization leads to 
a strong increase in NV'' signal for higher excitation powers. 
Ratio data points for which the NV" PL is not detectable 
have not been plotted. 

tocliromism, hypothesizing that an NV^ center can 
be discharged by high power excitation hght, re- 
sulting in an NV° center and a free electron. Fig. 
[7] shows the normalized PL intensity of NV" and 
NV~ with respect to excitation power, for a con- 
stant ion implantation density (F = 5 x 10^'^ ions 
cm^^). The normalized NV° PL intensity increases 
with increasing excitation power. At the powers used 
in our work (< 5 mW), PL saturation cannot be ob- 
served. For low excitation power, the NV" emission 
becomes negligible, as the vast majority of NV cen- 
ters is negatively charged, whereas with increasing 
excitation power, the fraction of neutral NV centers 
increases. A decrease of the signal strength of NV~ 
PL (normalized to laser power) occurs in the same 
way as NV*^ PL rises. Assuming that the number 
of NV centers is proportional to the photolumines- 
cence intensity multiplied by the radiative lifetime 
for each charge state (negative: 12 ns, neutral: 20 ns 
[T]), the total number of NV centers in the laser spot 
remains constant (within 4 %) for the range of ex- 
citation powers probed. This derivation is only ap- 
proximative, as it assumes the same excitation rate 
for both charge states, which is not necessarily the 
case at the given excitation wavelength. But it hints 
that conversion occurs only between the two charge 
states of the NV center, whereas optical conversion 
of the NV^ center into non-radiative centers, if oc- 
curing at all, is minor. With an excitation energy of 
2.33 eV of our system, this ionization process can 
only occur via a resonant two-photon process by sub- 
sequent excitation and ionization of the NV~ center 
(ground state 2.58 eV below conduction band [34]). 
The build-up time for the NV° PL signal has been 
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Fig. 8. PL ratio of NV~/NV'' versus implantation flucnce 
and laser excitation power (logarithmic scales). The figure 
shows clearly that a high ratio can be reached even for medi- 
um-high implantation densities, indicating temporal stability 
of the NV~ centers in high-density ensembles. Data points 
for which the NV'' PL is not detectable have not been plot- 
ted. In the figure, the surface has been obtained from the 
data points by spline interpolation. The black squares show 
the actual data points recorded. The PL linewidth (not plot- 
ted) remained constant at all laser powers within the accu- 
racy of our detector (0.1 nm). 

reported as on the order of hundreds of ^s [53] . At 
our peak excitation power, a conversion efficiency of 
up to 20 % could be achieved, ranging similarly as 
reported by Manson et al. [23j . The photoionization 
conversion efficiency is limited by spontaneous back- 
reaction of a conduction band electron and NV° cen- 
ter into an NV~ center, leading to an equilibrium 
for continuous-wave excitation. 

As the charge conversion process implies a change 
of resonance frequency and a loss of the electron 
to the conduction band with subsequent fast in- 
traband decoherence, QIP schemes involving single 
NV^ emitters would most likely fail in the case of 
light induced charge conversion. Schemes that in- 
volve ensembles of NV centers suffer a fidelity de- 
cline from a partial charge conversion of the con- 
tributing centers. It is therefore vital for QIP with 
NV~ ensembles to know the conditions for which 
the charge conversion of NV~ can be minimized. For 
the off-resonant Qmem scheme mentioned above, it 
is worth noting that the collective coupling strength 
of a high density NV^ ensemble allows for control 
laser powers much lower than the laser powers used 
here for excitation. Moreover, the two-photon charge 
transfer process is likely to be strongly suppressed 
for off-resonant frequencies. 

It is interesting to relate the photochromic behav- 
ior of the NV center to the ion implantation density. 
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Fig. 9. PL peak widths of NV°, NV~ and GRl with respect 
to irradiation fluence. An increase of the irradiation induced 
damage of the center environment leads to a stronger broad- 
ening of the luminescence line. 

Fig. [8] shows a systematic analysis of the NV~/NV° 
ratio as a function of implantation fluence and ex- 
citation power. As discussed above, at high implan- 
tation densities, the fraction of negatively charged 
NV centers decreases due to a lack of free electrons. 
Therefore, for very high implantation densities the 
NV° emission cannot be suppressed and the tempo- 
ral stability of the NV~ center is not guaranteed. 
On the other hand, for lower implantation densities, 
the presence of NV° can be completely suppressed 
within the accuracy of our optical detection appa- 
ratus (at a signal to noise ratio for NV~ PL of over 
100). Our systematic analysis (figure [8|) shows that 
such a regime can even be reached for medium ion 
implantation fluences of up to 5 x 10^^ ions cm~^. 
Therefore, the technique of ion implantation of ni- 
trogen rich HPHT diamond can be deployed to cre- 
ate high densities of NV~ centers without a de- 
tectable impact on their temporal stability. 

3.4. Inhomogeneous Broadening 

The linewidth of the NV^ PL peak in our sample 
is large, despite the cryogenic temperatures. This is 
typical for NV centers created by ion implantation 
in impurity-rich HPHT diamond. The peak widths 
measured (full width at half maximum, FWHM) 
monotonically increase from AA — 0.66 nm {/S.v ~ 
480 GHz) to AA = 2.7 nm (Ai/ ~ 2 THz) as a 
function of implantation fluence (see figure [9|) . We 
attribute this increase to the implantation induced 
damage to the crystal structure and an increase in 
stress, which is not fully repaired due to the com- 
paratively low annealing temperature. On the other 



hand, we could not detect a change of the NV~ 
peak width due to a variation of the excitation 
power. This is the case for all implantation regions, 
such that illumination induced broadening does 
not occur at any extent of implantation induced 
damage. We therefore conclude that the increase of 
free charge carriers due to high power laser illumi- 
nation (1.6 kW cm~^) has a small impact on the 
inhomogeneous linewidth compared to the broaden- 
ing mechanisms prevalent in high nitrogen density 
diamond. 

A large inhomogeneous broadening of the NV~ 
center ZPL has been reported for a great variety of 
samples [1]. Recent research has shown a great in- 
terest in a reduction of this broadening, as a large 
inhomogeneous linewidth of an optical emitter lim- 
its its usability for certain QIP schemes. For exam- 
ple, a quantum memory that stores a photon in the 
excited state of an ensemble of NV^ centers would 
dephase within picoseconds. Thus the controlled re- 
versible inhomogeneous broadening (CRIB) scheme 
[37|38| could not be realized in practice without ad- 
ditional efforts, e.g. spectral hole burning. Likewise, 
an FIT Qmem would require huge control pulse 
energies to achieve a usable transparency window 
(fully separated absorption peaks), which implies 
intensities above the material damage threshold. It 
appears that, with current characteristics of high 
density NV samples, resonant Qmem schemes can 
only be applied to single NV centers, where inhomo- 
geneous broadening is not present and spectral dif- 
fusion can be limited. For example, it has recently 
been shown [B] that a great reduction of the nitro- 
gen density by a careful choice of the raw sample 
material and deliberate implantation of single nitro- 
gen molecules can lead to NV" lincwidths of about 
Ai/ ^ 10 GHz. While this technique can create sin- 
gle NV centers, it is unclear whether high density 
NV ensembles with as little inhomogeneous broad- 
ening can ever be produced. It has to be reminded 
that on-resonant schemes as the FIT scheme require 
equally high optical densities {D as the here 

presented off- resonant scheme |33|39|40] . Moreover, 
a linewidth of 10 GHz is still two orders of magnitude 
larger than the natural linewidth (83 MHz [1]), such 
that for CRIB with spectral hole burning, the re- 
quired optical density equally increases by this large 
factor. 

The off-resonant Qmem scheme therefore appears 
advantageous, as inhomogeneous broadening of the 
excited state | 2) causes no difficulties. The ma- 
jor condition for off-resonant storage is a detuning 
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much larger than the photon bandwidth and the ex- 
cited state hnewidth. Its dcphasing time is subject 
to inhomogcneous broadening of the storage state 
only, not of the excited state. From the spectra 
recorded and using curve fitting, the inhomoge- 
neous linewidth can be estimated as Ai^ia ~ 2 THz, 
which corresponds to a local phonon dephasing 
time of 500 fs, which interestingly is a good order of 
magnitude shorter than the bulk phonon lifetime in 
diamond |41|42| . This would lead to rapid dephas- 
ing if this particular sample was used to build a 
PSB Qmem and the Qmem Q-factor would be low 
(Q = Aj/i3 6 ~ 1.5). However, no temporal analysis 
of the PSB lifetime has been published, such that 
the contributions of natural linewidth (leading to 
decoherence due to phonon decay) and inhomoge- 
neous broadening (leading to dephasing) are not 
clear. The latter could possibly be minimised for 
improved sample creation methods. As mentioned 
above, the PSB storage scheme is generic and can 
be applied to other color centers, including those 
with longer lifetimes of the local phonon mode. 

4. Conclusion 

We have analyzed the photoluminescence proper- 
ties of a type-lb diamond sample, which has been 
systematically He implanted over a wide range of 
fluences. At low and medium fluences, the PL signal 
from NV centers increases with the fluence, whereas 
at high fluences absorption and competing processes 
lead to a decline of the NV luminescence. The ra- 
tio of NV~/NV° fluorescence decreases for high im- 
plantation densities, as the density of residual nitro- 
gen atoms which act as electron donors decreases. 
We have analyzed this ratio with respect to irradi- 
ation density and excitation power. The concurrent 
examination by ion implantation fluence and exci- 
tation power reveals suitable conditions for which 
the NV~ center is a temporally stable absorber in 
the solid state, as the charge-conversion due to pho- 
toionization can be strongly suppressed. The results 
have been reviewed in the context of QIP. 

We have introduced a generic scheme to store 
broadband photons in ensembles of color centers in 
diamond via an off-resonant A-transition. A phonon 
sideband is used as the storage state, implying a 
short coherence time, which however is compensated 
by a large allowed bandwidth and ultrafast opera- 
tion. Due to inhomogcneous broadening, the only 
practical quantum memory type suitable for NV 



ensembles is the off-resonant (Raman) type. On- 
resonant schemes such as EIT arc more practical to 
implement with single NV centers, but require mi- 
crocavities. 
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